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We illustrate the contemporary application of Dyson-Schwinger equations using 
two examples: the calculation of pseudoscalar meson masses, an associated model- 
independent mass formula and the approach to the heavy-quark limit; and the 
study of nucleon observables, including a calculation of its mass, M, via a covariant 
Fadde'ev equation and an estimate of pion-loop contributions to M. 

1. Introductory Remarks. The Dyson-Schwinger equations (DSEs) in 
quantum field theory are an analogue of the Euler-Lagrange equations in clas- 
sical field theory. They are an enumerable infinity of coupled integral equations 
whose solutions are the n-point Schwinger functions (Euclidean Green func- 
tions). These Schwinger functions are also the matrix elements estimated in 
numerical simulations of lattice-QCD. The coupling between equations poses a 
challenge, of course: it necessitates a truncation in order to define a tractable 
problem. A weak coupling expansion provides one systematic method and 
reproduces perturbation theory. However, it also makes nonperturbative phe- 
nomena inaccessible, and something else is needed for the study of strongly 
interacting systems and bound state phenomena. 

This situation is familiar from many body theory where the Hartree-Fock 
truncation often yields reliable information. There is an analogue in QCD: 
the renormalisation-group-improved rainbow-ladder truncation, which is phe- 
nomenologically efficacious, as is clear from Ref. jjj wherein the power of a 
single-parameter model of the infrared behaviour of the effective quark-quark 
interaction is illustrated.^ The successes and failures of such a model can be 
understood once it is appreciated that the rainbow-ladder truncation is the 
leading order in a systematic, Ward-Takahashi identity preserving l/iV c -like 
expansionp But this demonstration is a single, small step toward a rigorous 
foundation for contemporary DSE modelling. 

One of the beauties of a model is that its simplicity makes possible a rapid 
comparison between theory and experiment. Another is that it can be wrong: 



a The ultraviolet behaviour is fixed and model-independent because this truncation of the 
DSEs reproduces perturbation theory. 
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an Ansatz is intuitively motivated and explored, and, if it is widely successful, 
its failures can point to unanticipated phenomena. Irrespective therefore of the 
difficulties that remain in providing a rigorous foundation for the application 
of DSEs in QCD, their phenomenological application plays an important part 
in attempts to elucidate nonperturbative phenomena. 

This is nowhere clearer than in the study of dynamical chiral symmetry 
breaking (DCSB). The DSEs provide the simplest and most direct means of 
understanding the dichotomy of the pion as both a Goldstone mode and a 
bound state of a massive dressed-quark and -antiquarkJj'Q and the result is 
model-independent. The analysis relies on the interplay between the QCD gap 
equation, the DSE for the dressed-quark propagator, and the inhomogeneous 
Bethe-Salpetcr equation (BSE) for the axial-vector vertex. This interplay is a 
consequence only of the axial-vectoi^Ward-Takahashi identity. Furthermore, a 
key quantitative conclusion follows QQ in order to reproduce observed charac- 
teristics of the spectrum, the kernel in the QCD gap equation must exhibit a 
significant enhancement on the domain Aq CD < k 2 < 2 GeV 2 . Identifying the 
origin of that enhancement; i.e., whether it is a feature of the dressed-gluon 
propagator alone or of the contraction of this propagator with the dressed- 
quark-gluon vertex, is currently an important focusQ a 

In the last decade the use of DSEs in QCD has attracted renewed interest 
and they have been applied to a broad range of phenomena, as is clear from the 
detailed summaries in Rcfs. ||,f[0]]. The approach has moved far beyond the 
calculation of indigent approximations to f„ and the vacuum quark condensate. 
Herein we illustrate this by focusing on two topics: pseudoscalar meson masses 
and their evolution with the current-quark mass; and a description of the 
nucleon, its interactions and the calculation of its mass. 

2. Pseudoscalar Meson Masses. Meson masses can be calculated by 
solving the renormalised homogeneous BSE: 

[T H {k; P)] tu = [ A [ XH (q; P)]„ K™(q, k; P) , (1) 

Jq 

where: r,. . . ,it represent colour-, Dirac- and flavour-matrix indices and H iden- 
tifies the meson under consideration; P is the total momentum and P 2 = —m 2 H 
is the eigenvalue condition for a solution; xn(q', P) '■= S(q+) r#(g; P) S(q-) is 
the Bethe-Salpeter wave function, with Tn{q]P) the fully-amputated Bcthc- 
Salpeter amplitude and S = diag(5' M , Sd, S s , . . .) the dressed-quark propaga- 
tor flavour matrix, q + = q + r)p P, g_ = q — (1 — rjp) P; [] and K(q, k; P) is 

b 'i~lP E [0, 1] is the momentum partitioning parameter. It appears because in a Poincare co- 
variant approach the definition of the relative momentum is arbitrary. Observable quantities 
must be independent of rjp. 
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the renormalised, fully-amputated quark-antiquark scattering kernel, which is 
two-particle-irreducible, with respect to the quark-antiquark pair of lines, and 
does not contain quark-antiquark to single gauge-boson annihilation diagrams, 
such as would describe the leptonic decay of a pseudoscalar meson. In Eq. (|j) 
J A := J A <i 4 g/(27r) 4 represents mnemonically a translationally-invariant reg- 
ularisation of the integral, with A the regulation mass-scale, and in this 
particular case the r.h.s. of the renormalised equation is cutoff-independent. 
(A — > oo is the final step in all calculations.) Equation (|l|) is obtained from 
the inhomogeneous BSE by equating meson pole residues. 

For a pseudoscalar meson the solution of Eq. (|l|) has the form 

r H (q;P)= (2) 

T% [iEh (q; P) + 7 • PFh (q;P)+-fqG H (q;P) + <J^q^ H H (q; P) 

where T H is the flavour matrix that specifies the mesonic channel under consid- 
eration; e.g., T K+ = (1/72) (A 4 + iA 5 ), with = 1, ... ,8} the Gell-Mann 
matrices. The amplitude is canonically normalised by requiring that the bound 
state contribution to the fully-amputated quark-antiquark scattering matrix: 
M = K + K{SS)K + ..., have unit residue. (See; e.g., Ref. |].) 

Using: the inhomogeneous BSEs for the axial-vector and pseudovector 
vertices; the dressed-quark DSE; and the fact that a nonperturbative Ward- 
Takahashi identity preserving truncation of the DSEs is possible, it was shown 
in Ref. || that, for flavour nonsinglet pseudoscalar mesons, 

f H m 2 H = M C H ri , (3) 

with M% := tr flavour [M (c) {T H , (T H )*}] , where M (0 = diag(m£,m^ m<, . . .) 
and (•)* indicates matrix transpose, so that M\j is proportional to the sum 
of the constituents' current-quark masses. This model-independent identity is 
valid for all current-quark masses, irrespective of their magnitude, and there- 
fore provides a single formula that unifies the light- and heavy-quark regimes. 

In Eq. (||), fu is the leptonic decay constant, which the derivation proves 
is given by 

f H P» = Z 2 J itr [(T")* 7 57„Xh(<z; P)\ , (4) 

where Z 2 = Z 2 (£, A) is the dressed-quark wave function renormalisation con- 
stant, with £ the renormalisation point. This multiplicative factor of Z 2 on the 
r.h.s. ensures that £jj is gauge-invariant, and independent of £ and A; i.e., that 
it is an observableo Equation (0) yields fn as the pseudovector projection of 
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the meson's Bethe-Salpeter wave function at the origin in configuration space; 
i.e., this equation provides a field theoretical analogue of the "wave function at 
the origin," which describes the decay of bound states in quantum mechanics. 
The remaining term in Eq. (||) is 
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(T H f ^XH(q;P) , (5) 



where Z4 = Z±{£, A) is the dressed-quark mass renormalisation constant. The 
gauge dependence of Z4 is precisely that necessary to ensure that the r.h.s. 
of Eq. (||) is gauge invariant; its cutoff dependence ensures that the r.h.s. is 
independent of the cutoff; and its renormalisation point dependence ensures 
that the product on the r.h.s. is independent of the renormalisation pointll r H 
is the pseudoscalar projection of the meson's Bethe-Salpeter wave function at 
the origin in configuration space. 

In asymptotically free theories the chiral limit is unambiguously definedtHa 
by setting rh = 0, where ah is the renormalisation point independent current- 
quark mass. In this limitcl 

r k := lim n t h = -7— > ( 6 ) 

where fjj a is obtained by taking the chiral limit in Eq. (|J) and 

,A 

-(qq)° c = Z A N c / tY D [S A=0 (q)} , (7) 



with Sm=o obtained as the chiral limit solution of the dressed-quark DSE: this 
is the gauge-invariant and cutoff-independent expression for the vacuum quark 
condensate. Using Eq. (^), Eq. (||) yields 

fk m% = - (m% + m C h ) (qq) ° + O (m 2 fl j2 ) , (8) 

with fi t 2 labelling the flavour of the dressed-quark constituents; i.e., as a corol- 
lary, Eq. (||) yields the so called Gell-Mann-Oakes-Renner relation. 

As remarked above, Eq. (^) is also valid for arbitrarily large current-quark 
masses and an analysis of its heavy-quark limit is facilitated by writing 
H111 — (Mq+Eh), where Mq is a constituent-heavy-quark massE3 and Eh is 
a "binding energy." Following this the dressed-propagator for the heavy-quark 
constituent takes the form 




Sq(k + P) = L J——L— + O [ ) , (9) 
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where k is the momentum of the lighter constituent, and the canonically nor- 
malised Bethe-Salpeter amplitude can be expressed as 



T H {k;P) = V^Fr<°°(fc;P), (10) 

where r^°°(fc; P) is pxiintwise finite in the limit m H — > oo. Using Eqs. (||) and 
© in Eq. (§ yieldsU 

/h = -^=, (11) 

with c H a calculable and finite constant, which reproduces a well-known con- 
sequence of heavy-quark symmetry. Applying the same analysis to Eq. (||) one 
fmdsta 

r£ = ^ (12) 

and this, along with Eq. ([ll]) in Eq. (||), provesEBIll that in the heavy-quark 
limit 

mH = ^fM C H ; (13) 

i.e., that the mass of a heavy pseudoscalar meson rises linearly with the mass 
of its heaviest constituent. 

It has been showr£3 that Eq. ( pi] ) is not valid until current-quark masses 
m > mi,. The c-quark lies well-outside this domain; e.g., if a constant of 
proportionality is chosen so as to reproduce the value of then fjj obtained 
from this formula is ~ 40% too large. This is consistent with the calculated 
magnitude of the violations of heavy-quark symmetry in b — > c transitions 
(< 30%). 

Using the results reported in Ref. El , obtained using the renormalisation- 
group- improved rainbow-ladder truncation of Ref. H] , one can study the evo- 
lution of meson masses as the current-quark mass is increased. Consider first 
pseudoscalar mesons whose constituents have equal current-masses, for which 
the calculated evolution is described by the interpolating formulaEl 

mjj, =pVx + ~fX, (14) 

where = 1.04 GeV, 7 = 0.21 GeV and X = m f /A Q CD, with C = 19GeV and 
Aqcd = 0.234. This formula was determined via an unconstrained fit to the 
masses calculated by solving the Bethe-Salpeter equation. A comparison with 
Eq. ((§) shows that := /3 2 /(2Aqcd) = — (??)°//h anc l> usm g the model's 
calculated value of fu = 0.088 GeV, one infers a value of (qq)® = (—0.26 GeV) 3 
from this correspondence. That can be compared with the value (— 0.27) 3 
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calculated directly in this model from Eq. (Q). This near equality indicates 
that the interpolating formula in Eq. ( |l4| ) can provide reliable estimates. 

Equation ( |l4| ) indicates that a flavour nonsinglet, fx = f = pseudoscalar 
meson with a mass m# = 1 GeV would be composed of quarks with mass 
rhf = 0.32 Aqcd = 2.3 rh s . At this current-quark mass, which corresponds to 
X = 0.68, the \f~X~ term still provides 86% of the meson's mass. Thus one 
remains well away from the linear trajectory, in spite of the fact that in the 
neighbourhood of X = 0.68 a tangent to the curve in the left panel of Fig. 2, 
Ref. jl| , is nearly indistinguishable from the curve itself within the resolution 
of that figure. Furthermore one finds easily from Eq. (IIJ) that 



mf=2ms ^ 22 



in agreement with the result obtained in recent lattice simulations^ Hence, 
in addition to being phenomenologically efficacious, the renormalisation-group- 
improved rainbow-ladder truncation of Ref. |L4| predicts a mass-evolution that 
is confirmed by lattice simulations. These results support a scenario of DCSB 
in which the vacuum quark condensate is large; i.e., _Bq GcV ^> fjj a - 

Reference [Q also provides an interpolation of the kaon-like u-q trajectory: 

m H K = 0.083 + 0.5 y/X + 0.31 X, (16) 

and while the BSE studies reviewed in Ref. have not directly addressed 
heavy-light nor heavy-heavy bound states (systems with rhf > 3.5 rh s have 
not been studied) one may, as a preliminary step, ask whether Eq. ) can be 



one 

V r^. 



used to obtain reliable mass-estimates via extrapolation? Using Eq. (lq) 
reproduces m D ~ 1.9 GeV, m B ^ 5.3 GeV with mJ GoV ~ 1.1 GeV, mjp - 
4.2 GeV, fund since these current-quark masses are in agreement with other 
estimatesta then Eq. ( |l6| ) can be a useful tool. In this application one finds 
that the linear term provides 50% of mo and 67% of tub- Thus, like the 
conclusion drawn on the validity of Eq. ( pT| ) for fjj, the heavy-quark limit in 
Eq. (|l^) is not valid until the current-quark masses satisfy m > nib, and the 
c-quark mass is much smaller than this lower bound. 

3. A Model of the Nucleon. The succesSpOf the rainbow-ladder trunca- 
tion in describing meson observables motivatesLj a treatment of the nucleon 
as a bound state of a dressed-quark and nonpointlike diquark via a covariant 
Fadde'ev equation. The feasibility of this approach was demonstrated in Ref. 
fill and the most extensive study to date is described in Ref. [|9|. The ap- 
proach assumes only that the colour-3 quark-quark scattering matrix can be 
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approximated by a sum of diquark pseudoparticle terms: scalar + pseudovec- 
tor + . . . , whose properties can be determined independently. The Fadde'ev 
equation then describes the nucleon as a quark-diquark composite, which is 
bound by the repeated exchange of roles between the dormant and diquark- 
participant quarks, and the complete nucleon amplitude is a sum of three 
terms: 

# = *l + #2 + *3, (17) 

where the subscript identifies the dormant quark and; e.g., \Pi,2 are obtained 
from \J/ 3 via a cyclic permutation of the indices: *i = ^(l^li), ^2 = ^3(1312)- 
The simplest such model retains only the contribution of the scalar diquark 
to the quark-quark scattering matrix, in which case 

^ 3 {pi;ai,Ti;a,T) = 
e ClC2C3 A 0+ (K) [T 0+ (^p [12] ;K)Z\Z ^(/;P)u(P) , (18) 

where: (ry ■ P + M)u(p) — 0, with P = pi + P2 + Ps =: P{\23] the nu- 
cleon's total momentum and M its mass; e clC2C3 is the colour-singlet factor; 
K = pi + p 2 =: P{i2j, P[i2] ■= Pi ~ P2, £ = (P{i2\ - 2p 3 )/3; (ai,Ti) are 
the quark spinor and isospin labels, and (a,r) are those of the nucleon. In 
Eq. (|l8|), A 0+ is the pseudoparticle propagator for a scalar diquark formed 
from quarks 1 and 2, and T + is a Bethe-Salpctcr-likc amplitude describing 
their relative momentum correlation. Both these quantities are determined 
by studying the quark-quark scattering matrix. The remaining element, 1P3, 
is a Bethe-Salpeter-like amplitude that describes the relative momentum cor- 
relation between the third quark and the diquark's centre-of-momentum. It 
satisfies a renormalised Fadde'ev equation, which in the isospin symmetric limit 
assumes the form 

^{k;P)u(P) = -2 / A 0+ (^)r 0+ (fc + ^/2;X) 
Jt 

x S(4x)* f 0+ (t + k/2; -K k ) S{h) P) u(P) , (19) 

with K t = -£+(2/3)P, £ cx = -i-k-P/3, l x = £ + P/3. The general solution 
for a positive energy nucleon takes the form 

Mb p) = h{i;P)t- ^ (ii- 1- e- Pi) W\ p) , (20) 

with P 2 = —1 and where, in the nucleon rest frame, fi t2 describe, respectively, 
the upper/lower component of the dressed-nucleon spinor. 
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([18]) has been used successfulb/-.tQ 

rs.BB 



The nucleon amplitude in Eqs. (|l7j), ( p"8| 
calculate a wide range of lcptonic and nonlcptonic nucleon form factor 
In those calculations, solving the Fadde'ev equation was side-stepped by em- 
ploying simple parametrisations for the functions in Eqs. (Eq), (^o|): 



A 0+ (K 2 ) 



1 



'o+ 



T 0+ (k;K) = Ci lb iT 2 T{k 2 /uul + ). 
f 1 (£;P) = ^ r f(l 2 /^), 



(21) 
(22) 
(23) 



with / 2 = r/i , where R is a constant of proportionality that gauges the relative 
importance of the lower component of the nucleon spinor. In these equations 
C = 7274 is the charge conjugation matrix, T(y) = (1 — e~ v )/y, and Af^, A/"o+ 
are calculated, canonical normalisation constants. 

The parameters in this model are u>q+ , m + , u),p : d + =1/ loq+ measures 
the quark-quark separation in the diquark; £ + = l/mo+ is the diquark corre- 
lation length or mean free path; and = measures the quark-diquark 
separation. | Their values have been determined by requiring a good impulse- 
approximation fit to the proton's charge form factor on Q 2 E [0, 3] GeV 2 and 
this procedure yields 
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0.64 , 
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= 1 
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0.29 


0.31 
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= 0.65 


0.21 


0.72 


0.63 
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= 0.65 


0.95 


0.28 


0.32 



(24) 



The scalar diquark parameter values determined in these unconstrained fits 
are within 10% of those obtained in the BSE studies of Ref. 

Using the values in Eq. ( p4] ) a wide range of observables can be calculated 
and herein we exemplify the results via 
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Calc. R=1 
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(0.81) 2 


-(0.34) 2 
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2.82 
2.85 


-1.62 
-1.63 


-0.57 
-0.57 


14.6 
14.5 


1.27 
1.12 



(25) 



and the neutron electric form factor in Fig. One significant feature, apparent 
in Eq. (p5|), is that \^ n /^p\ > 0.5, which is only possible because the impulse 

C A description of the nucleon in this form can only be internally consistent if d + < and 
£q+ < d^p; i.e., the diquark is smaller than the nucleon and can't propagate over distances 
larger than the nucleon. 
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q 2 (GeV 2 ) 



Figure 1: Ne utron electric form factor calculated with thip* different values of R . Data 
from Ref. |23| , extracted using the Argonne V18 potentialcj is very sensitive to the 
strength of /2. (NB. The contribution of the pseudovector diquark was not included in 
these calculations.) 

approximation explicitly includes diquark breakup contributions. Another is a 
prediction for the ratio ^, p G p E {q 2 ) / G p m (j£) that is in semi-quantitative agree- 
ment with recent results from TJNAFtJ As shown explicitly in Ref. ||, the 
impulse approximation requires five terms when the diquark is composite, a 
fact also appreciated in Ref. p^j. It also requires an explicit form for the 
dressed-quark propagator but that is well-known from studies of meson ob- 
servables and is given; e.g., in Refs. ]2C| , |2l| . 

Hitherto the choice of R is arbitrary. However, its value is fixed by solv- 
ing the Fadde'ev equation. To exemplify that we have solved the Fadde'ev 
equation, Eq. (|l9|), with the values of w +, rn + in Eq. ( |24| ) and assuming that 
f i (£; P) — f i {jP). This yields the results in Table [l| As anticipated in Ref. 
pl[ , a reduction of ~ 30 - 40 % is required in the calculated mass for agreement 
with experiment and here the inclusion of a pseudovector diquark can helpE2 

To explore that we added such a correlation to Eq. (Q) : 

^ 3 1+ =e ClC2C3 All(K)t ltl Cer 1+ ( 1 - P[12] )Al(e;P)u(P), (26) 
where: {t i= ( +! o,-) = (r + t 3 )/V2, n, (r - t 3 )/V2}, r = diag(l, 1), are the 
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Table 1: Nucleoli mass and R = /2//1 ratio determined by solving the Fadde'ev equation. 

, ^iij 2 are the widths of a least-squares fit to /i, fa assuming they are pointwise well- 
approximated by !F(£ 2 /uj^ f ). [They actually fall faster with increasing k 2 but these widths 
are nevertheless a useful guide for comparison with Eq. (pj|).] Rows 1 and 5 give the results 
obtained with only a scalar diquark while the others were obtained with the inclusion of 
a pseudovector diquark correlation. All dimensioned quantities are given in GeV. In these 
calculations the dressed-quark's Euclidean constituent-q«<jtrk mass is = 0.33 GeV, which 
was fixed in independent studies of meson observablesE3 Using the parameters in the last 
row, an analogous Fadde'ev equation for the A(1232) yields a mass = 1.23 GeV. 
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0.64 


0.40 


0.82 


1.11 


0.29 


0.31 


1.14 


0+ & 1+ 
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0.64 






2.11 


0.34 
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0+ & 1+ 
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0.64 
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1.26 


0.31 


0.32 


1.20 


0+ & 1+ 


0.42 


0.64 
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0.86 


1.68 


0.34 


0.33 


1.27 


0+ & 1+ 


0.42 


0.64 


1.09 


0.86 


0.46 
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0.35 
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symmetric isospin-triplet matrices, 

Al(£; P) = a\(£ 2 ) + a 2 (£ 2 ) 7s7 • Vl v , P = 1 , (27) 
with a\ = a\ = a\ , a\ = a\ = a 2 in the isospin symmetric limit; and 

{^ + ^) -^HK 2 /^ + ), (28) 

r 1+ (fc 2 ) = -^-^(fc 2 /^ + ); (29) 

and extended Eq. ([l9]) to include the coupling to this channel. 

The results obtained in this case are also presented in Table [l| In the 
upper panel we chose the value of mi+ /m + from Ref. [^2| and, initially, uj\+ — 
lu + , to find that the pseudovector diquark provides additional attraction and 
reduces the calculated mass by 11%. Reducing W]+ increases the pseudovector 
diquark coupling via an increase in gi+ := hence the calculated mass is 

reduced: a 40% reduction in loi+ reduces M by 10%. Reducing ui Q + by the same 
amount has very little effect, reducing M by only an extra 2%. However, in 
agreement with intuition, decreasing the diquark masses reduces the calculated 
nucleon mass: a 21% decrease yields a 19% reduction in M. In the lower panel, 
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we expanded the system of coupled equations to include an analogous Fadde'ev 
equation for the A. The exploration of this system made clear that an increase 
in g x + via a decrease in uj x + is not the only way to reduce the nucleon mass. 
Increasing uj x + also increases the support of the integrand that describes the 
pseudovector's binding contribution and, for u>i+ greater-than a certain value, 
that increase more than compensates for the concomitant decrease in g x + , 
yielding increased binding and a lower nucleon mass. The table makes amply 
clear that an internally consistent description of the nucleon and A is possible 
using only scalar and pseudovector diquark correlations, just as found in Ref. 
p9| . Furthermore, it is clear from the table that a well-constrained scalar 
diquark model should employ R in the range ~ 0.5 — 0.7. 

A question that remains unaddressed is the role of pion loops. The on- 
shell ttNN coupling is strong and hence it is conceivable that such loops might 
generate large self energy corrections to the nucleon's mass. We have made an 
estimate using a model DSE for the nucleon self energy: 

G- 1 (p) = i 1 -p + M + Z J ^(p~k)gl NN {{ P -k) 2 ,k 2 ) l5 G{k) l5 , (30) 

where A^(£) = l/(£ 2 +m 2 ) and g-nNN{{p—k) 2 , k 2 ) is the momentum-dependent 
ttNN coupling. The t = — (p— fc) 2 -dependence of this coupling was calculated 
for on-shell nucleons in Ref. j2^] , with the result 

g*NN(t, k 2 = —M 2 ) a - gff "* , A, = 0.96 GeV, (31) 

(i - t/Aiy 



where g-^NN is given in Eq. (|25|). This is not quite sufficient for our present 
purpose because the nucleon in the loop is not on-shell. Therefore to complete 
an estimate we employ a simple product Ansatz: 

g,NN(p; k) = (i - |p2 - fc2 | /A 2 )2 (1 + p + M 2| + \ k 2 + m 2 \)/AI) 2 ' (32) 

to approximate the angular-average of the coupling that is active in the integral 
equation. With the parameters in the last row of the upper panel in Table [l| 
we find that the pion loop adds lOMeV to the nucleon's mass; i.e., it provides 
only a 1% increase. (The detailed form of Eq. (|3l]) is not important but the off- 
shell suppression is.) In addition, if we define M m7r to be the loop-corrected 
nucleon mass obtained using a particular value of m n , then M m7r decreases 
with increasing m ff , in qualitative agreement with chiral perturbation theory. 
(See; e.g., Ref. jj.) 

d The positive sign before the integral on the r.h.s. is correct and entails that both the vector 
and scalar self energies are positive. 
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Epilogue. Two short summaries are all we have room for here. In large part, 
the light-quark meson sector is well understood. The renormalisation-group- 
improved rainbow-ladder truncation provides reliable information in many 
channels and where it doesn't the reasons why are understoodH Understanding 
the nucleon: it's mass and interactions, is a contemporary focus and progress is 
rapid. Success there will open the way for many new phenomenological applica- 
tions, such as the calculation of quark distribution functions, which are merely 
parametrised in analyses of deep inelastic scattering. That the DSEs can pro- 
vide valuable insight here_is demonstrated by a calculation of the valence quark 
distribution in the pionO which itself is measurable given a high-luminosity 
electron-proton collidcrEJ Analogous, detailed Fadde'ev-equation-based stud- 
ies of the A are also beginning, as illustrated herein and in Ref . |3l| . 
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